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T.HE  QUASIGEOID  A»  THE  STA1DAHD  ELEVATIONS 
-  EAST  GEBiANY  - 

^FollcnTing  is  a  translation  of  an  article  by  K.  Mueller, 
German  Chamber  of  Technology,  scientific  contributor 
of  the  Leipzig  Geodetic  Service,  in  the  German  language 
periodical  Vei''messungstechnik  (^odesy).  East  Berlin, 

No.  3,  Mar cK T960 , " ^age s '  6T-o^^7 

(Conclusion  from  No.  2,  196C) 


6.  Brims'  Theorem 

faking  point  A  on  geoid  W  =  then,  according  to  equation 
(29a),  the  following  equation  holds  for  point  B  of  the  topographic 
earth  surface  B 


W3  ~  W  -  ~  j  gdh, 
A 


(31) 


For  the  same  noint  B,  the  normal  potential  U_  then  looks  as  fol- 

lOV^S!  ‘  B 


Ug  ~  -  —  V  V  dh  I 


f  r  dh  —  F  dH. 

0  o  B> 

taking  model  earth  as  reference  ellipsoid  then  is  (see  Fig.  3) 

B' 

According  to  (I6)  we  have  —  ^  dh  =  Ug,  — 

and  according  to  (13);  we  have  Up,  =  Wp  — 

From  (16)  and  (13)  follows  dh. ' 

3'  B 

Using  (31),  we  get  —  \'Y  dh  -  --  j  gdh 

o  '  ^ 

and  from  (32)  follows  g 


(32) 


h 


A 

v*  d  H. 


(33) 


B'. 


For  the  last  integral,  hoirever,  one  can  take  ^  as  constant  and, 
because  of  the  meaning  of  ,  v/e  can  set  up  the  following:  : 

B  -  B  '  . 

5  r  d  H  =  T  (  d  H  =  Y'  o 
B'  B'  '  .  ■  E  ,  , 

Thus  we  get  Ug  —  U  -  —  ^  gdh  — , 

°  k 


-  1 
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From  (31)  and  (3U)  we  get  the  following  by  forming  the  difference 


Wb  - 


u. 
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=  r 


(35) 


or 


f- 
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If  we  define  the  perturbation  function  Tb  as  follows: 


Tb  =  Wb  -  Ub 


(36) 

(37) 


and  if  we  consider  that  the  model  earth,  which  defines  the  quasi- 
geoid,  has  been  chosen  because  of  its  property  5  in  such  a  way 
that  TIq  -  Uq,  then  we  finally  get  the  well-knoim  theorem  by  Bruns, 
after  dropping  the  index  B,  using  ^  ^ 

,(38) 

But  here  ^  is  the  normal  gravity  on  the  physical  earth  surface 
and  ^  the  altitude  of  the  quasigeoid. 


7 .  Development  of  an  Integral  Equation  for  the  Altitude  of  the 
Quasigeoid  __  "" 

According  to  a  calculation  by  Levallois  l2/f  we  novi  want  to 
develop  an  integral  equation  for  ;  this  equation  was  obtained 
by  Molodjenski  /S7in  a  different  way  and  was  mentioned  by  Sagre- 

bin  /Io7. 

Let  us  start  from  Green's  theorem 

y  ill  —  V  j  y  ds  =  /Y  .  a  V  —  V.il  1/  d  r  ,  (39) 

^  iJ  ^  TTH 

where  I  is  a  volume  v>rith  the  surface^S,  and  Y  and  V  are  two 
harmonic  functions  of  space  and  ~  -f- 

is  the  Laplace  operator,  ^  "  ciZ. 

Here  we  take  for  Y  the  gravitation  potential  W  of  the  earth 
and  put  y  -  l;r,  where  r  is  the  distance  of  the  moving  point  P  on 
the  surface  S  from  a  fixed  point  M  on  this  surface.  Let  the  volume 
extend  into  infinity.  The  direction  of  the  outer  normal  n  of  the 
volume  "V  points  from  the  point  P  to  the  earth.  If  we  separate 
point  K,  which  lies  on  S,  from  the  points  of  the  volume  ’7^'  ,  by 

means  of  a  small  hemisphere  ^  ,  placed  on  S,  around  M  as  center 

with  infinitesimal  radius  S  f  then  the  inner  boundary  of  volume 
consists  of  the  surface  S,  let  this  surface  be  S'  —  minus  the  small 
circle  that  has  been  drawn  with  the  radius  S  around  M  —  and  of  the 
surface  of  the  hemisphere  around  M. 

In  this  region  "77  ,  W  and  l:r  are  harmonic  functions^  for 

them,  the  Laplace  expression  thus  becomes  W  -  0  and  Ti  -p-  -  0  . 
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The  right  side  of  Green's  theorem  (39)  disappears  too  and  one 
obtains  _ 


-i-  iL^  f* 

3 


o  I 

>.i  t. 


J-  lii  J  £  -  c 


For  the  surface  element  d  ^  one  obtains  in  polar  coordinates; 

d  ^  ^  sin  9  d  >  d/s  (1 

and  for  the  hemisphere  we  furthermore  get: 


-  a  4r- 


^  -h  iir,  , 


The  following  is  valid  for  and  A  of  the  hemisphere  A 
—  )i  t  y'  ='  /  ]i.  and  0  ^  A  ^  -A  . 

From  this  then  follovfs 
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For  lim  5  -  0  vre  get  from  (iiO)  and  (Uii) 
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^  ^  4“o/n  aJ  -  (2i, 

‘J  V'-'O  ^  ^ 

"r  ® 

Herein  noW  lim  W  -  Vf|,j  and  the  second  integral  becomes 

T<  ^ 7\ 


(h^) 


>.  r*  *'•  h 
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The  last  integral  in  (1^5)  vanishes  ’with  lim  C 
vfe  obtain  from  (ij5) 


M  Mi  .  I 

/  ,  (  w  r  —  1  a  W  ) 

^  ^  (  c>  n  .  r  3  n  ) 

Finally  one  obtains  from  this 

n  (  1 


dS 


-  0.  Using  this, 
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2rt  =  0. 
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11 
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(  £_  “  1  Q  W  ) 

I  (W  d  n  r  ’’TJ^Ia  )  dS, 

2~rr  U  (  ) 
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(h7) 


Applying  the  same  formula  to  the  normal  potential  U  —  this  is  per¬ 
missible,  since  U  satisfies  all  conditions  of  Green's  theorem  — 
one  then  obtains  in  an  analogous  way 

f  C'  ^  a  i  ) 

%  - _ L  \  )  ^  JL.  ~~  ^  cJ  U  )  <^3. 

^■*'1  <)  U  (u  a  n  r  “‘i]  n  )  (u8) 

S 


Comparing  the  formulas  (li?)  and  (iiB),  we  obtain 


-  h  - 


%  % 
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(W  — U)  r  —  1  (J  YJ  U)  j 

An  r  (a  n  d  n)  i 


(1^9) 


According  to  Bruns'  theorem  (38)  we  now  have  T  =  ^  -  Ya-U  and 

one  obtains  . .  ^ 

ot?  1 


\  u  Y  - 


n  pd  1 

1  \  V  L-  ^  di-I ) 

“ti  J  ,lf(^dn  r  (a  n  a  n  ) 


dS. 


(50) 


In  the  integrand  ^  \l/  cJ  n  and  £)  U/^  n  are  the  derivatives  of 
W  and  U  in  the  direction  of  the  normal  on  the  topographic  earth  sur 
face.  In  the  moving  point  P  of  the  physical  earth  surface  th^ 
latter's  normal  forms ^  with  the  directions  of  the  fields  g  or  il 
the  angles  ry  ort>,' '—cX  ,  respectively.  Then^  because  of  (i,)  or 
(13)  and  considering  the  relations  in  Figure  6,  we  get 


d  w 

-  0  y; 

cl  n 

“  a  H 

z 

d  n 

c,  H 

3  K  ;  __ 


d  n 


gp 


CCS  x 


d  V!  ~  _^U  ,  ^  H  =■•  —  Y'  p  COS  ex'  •“  —  f  p  COS  >X. 

Using  (51),  equation  (50)  then  goes  over  inoo 


(51) 


$  H  V'  =  1  )  1 

271  J  U 


b  1  , 


dS. 


(52) 

According  to  free_air  reduction  (lU),  the  following  equation  holds 
for  the  altitude  H  =  ^  of  P  over  the  model  earth  as  reference 
surfaces 


vp  =  \  q  -  2  ii-i 


H 


One  therefore  obtains 

(i 


) 

N 

!  ) 


(53) 


v  =  j^  \  ) 

2  71  U  f 
S 


V  V  r 


nj  n 

^  \  [ep  -  Vi,  U  -  2  ii 


>  COS  cx 


dS, 


(5ii) 
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If  we  want  to  introduce  into  this  equation  the  "free  air  anomaly" 

•4  g)  we  must  reduce  gp  with  the  altitude  of  P  to  the  quasigeoid, 
or  —  which  amounts  to  the  same  thing,  except  for  second  order  mag¬ 
nitudes  —  relate^  to  the  level  of  gp.  Using  (lli)  we  get 

4  g  =  gp  --  '^P  =  gp  “  ^  “I"  )  • 

Finally  we  obtain  the  desired  integral  equation 
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(56) 


This  equation  gives  the  quasigeoid  as  function  of  the  free 
air  anomaly  on  the  topographic  earth  surface  and  determines  the 
value  5  distance  of  the  quasigeoid  from  the  reference  ellipsoid. 

The  solution  of  the  integral  equation  (56)  should  be  possible 
by  iteration  if  one  knows  enough  values  of  the  gravitation  anomaly 
on  S.  A  solution  in  first  aporoximation  for  the  integral  equation 
(56)  one  obtains,  as  Molodjenski  shovred  in  his  work,  from  the  Stockes 


formula 

C 

0  M  =  R 

\  '\  (g  5'  (>• )  d’/'  .  -  y<: 

rnio 

J  J 

O  0 
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F(U') 


we  have 

-  a  mean  earth  radius; 

-  a  mean  normal  gravity  on  the  earth; 

-  the  measured  value  of  the  gravity  in  uhe  moving  point  P; 
the  normal  value  of  the  gravity  which  corresponds  to  the 
measured  value  g;  it  is  calculated  with  the  normal  alti¬ 
tude  H  and  the  latitude  ('*  of  the  point  P; 

the  angular  distance  of  the  point  P  from  the  fixed  point 
M  which  is  to  be  investigated;  it  runs  from  0  to  2  Tf  I 

-  the  azimuth  of  the  arc  MP;  it  runs  from  0  to  Ti'  ; 

-  the  yrell-knoTiTn  Stockes  function 


F  (W  )  =  1  sin  '>ij 

"T* 

—  3  cos'ji  In 


cosec 


T'  —  6  sin 

"T* 


2L  ^  1 


—  5  cos  'jj 


f 

^sin 


I 


sin 
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A  derivation  of  the  Stockes  formula  without  the  us0_of  spherical 
harmonics  (functions)  can  be  found  in  Baeschlin  ^12/ a  For  prac¬ 
tical  aoplications ,  this  formula  was  prepared  by  W.  F.  Jeremejew 

'  ... 

The  numerical  evaluation  of  the  integral  xs  made,  using 

the  same  process  as  P.  S.  Sakatow  uses  it  in  his  book  /ill/ 
the  calculation  ox  the  components  of  the  plumb  deviations.  Bet¬ 
ter  values  are  obtained  by  using  the  palettes  of  Jeremejew  and 
Jurkina  /iS/  ’i/ifith  which  one  can  calculate  X  up  to  a  few  dm. 
This  is  a  quite  satisfactory  result  for  the  geoidal  part  of  the 


altitude. 

We  must  stress  again  that  both  altitude  portions  can  be 
ca3.culated  exactly.  This  cannot  be  done  for  the  orthometric  alti¬ 
tudes. 


8.  Conclusion 

It  is  quite  justified  to  ask  for  an  exact  solution  for  the 
calculation  of  the  distance  of  a  point  on  the  physical  earth  sur¬ 
face  from  the  reference  ellipsoid.  One  must  consider  that  all 
measurements  of  geodetic  quantities  are  performed  on  the  physical 
earth  surface.  This  surface^  however,  is  not  suitable  for  relating 
on  it  the  measured  values  as  far  as  their  interdependence  is  con¬ 
cerned.  For  this  purpose,  the  measured  elements  of  geodetic  grids, 
e.g.,  base  lines,  astronomic  coordinates,  horizontal  directions, 
etc.,  have  until  nov^  been  reduced  to  sea  level,  i.e.,  to  the  geoid . 
From  the, geoid  the  quantities  vfcre  transferred  to  the  reference  el¬ 
lipsoid.  Distances  in  this  procedure  were  transferred,  under  con¬ 
servation  of  the  magnitude,  from  the  geoid  onto  the , reference  el¬ 
lipsoid,  neglecting  the  difference  in  curvature  of  both  surfaces. 
This  procedure  is  not  exact  mathematically  and  adds  further  uncer¬ 
tainties  to  the  uncertainties  connected  with  the  geoid,  sc  that 
thc3  value  of  precise  lueasureraents  can  be  strongly  reduced.  This 
procedure  therefore  must  be  rejected.  A  method  that  avoids  these 
errors  is  the  '’projection  procedure”  developed  in  the  USSR.  In 
this  procedure,  the  measured  geodetic  quantities  are  projected 
along  the  lines  of  force  to  the  quasigeoid;  they  are  then  pro¬ 
jected  along  plumb  directions  to  the  reference  ellipsoid.  As  Pi- 
zatti  suggested,  the  projection  is  tnus  perfox'med  in  two  steps. 

A  confrontation  of  both,  methods  is  given  by  hilan  Bursa '/l^/.  The 
Pizetti  projection,  with  the  quasigeoid  as  intermediate  sui-face, 
can  be  performed  without  the  use  of  any  hypothesis.  Here,  all  the 
quantities  that  are  normally  neglected  are  perfectly  knovm.  This 
clearly  shows  the  great  importance  of  the  idea  of  M.  S.  Molodjonski 

If  one  takes  into  consideration  finally  tnat  the'  quasigeoid 
coincides  on  the  free  oceans  Virith  the  geoid  of  -Listing  and  that  it 
touches  it  on  the  shores  of  the  continents  ^ —  i.,e . ,  here  only  the 
radii  of  curvature  differ  —  and  if  one  further  considers  that  the 
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quasigeoid  has  only  a  small  distance  from  the_geoid,  then  one  can 
understand  the  suggestion  made  by  Levallois  (_  9_7  "that  the  quasi¬ 
geoid  be  defined  as  earth  figure.  One  can  do  this  especially  since 
its  surface  can  be  determined  with  satisfactory  exactness  from  the 
integral  equation  that  was  developed  above,  if  one  only  knows  suf¬ 
ficient  values  of  the  anomaly  spaced  equally  over  the  vAole  earth. 
At  least  this  definition  would  be  no  less  conventional  than  the 
geoid  used  by  Listing,  and  the  equations  used  foi’  the  quasigeoid 
are  exact.  If  in  the  future  the  dimensions  of  the  model  earth 
(Kras sows ki  ellipsoid)  and,  with  this,  the  dimensions  of  the  quasi¬ 
geoid  should  change  slightly,  then  all  obtained  values  could  be 
transformed  exactly  into  the  new  values.  This  is  a  further  adv/'anb- 
age  of  the  quasigeoid. 

All  these  facts  caused  the  Soviet  geodesists  to  introduce 
for  the  USSR  the  normal  altitudes  and  thus  the  quasigeoid. 
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